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Abstract 

We study QCD in 1+1 dimensions in the large limit using light-front Hamil- 
tonian perturbation theory in the \/Nc expansion. We use this formalism to exactly 
compute hadronic transition matrix elements for arbitrary currents at leading order 
in l/Nc- We compute the semileptonic differential decay rate of a heavy meson, 
dr/dx, and its moments, M^v, using the hadronic matrix elements obtained previ- 
ously. We put some emphasis in trying to understand parity invariance. We also 
study with special care the kinematic region where the operator product expansion 
~1 — x~l)or non-local effective field theories ~ 1 — x ~ Aqcd/itiq) 

can be applied. We then compare with the results obtained using an effective field 
theory approach based on perturbative factorization, with the focus to better un- 
derstand quark-hadron duality. At the end of the day, using effective field theories, 
we have been able to obtain expressions for the moments with relative accuracy of 
0{A'QQjj/mQ) in the kinematic region where the operator product expansion can be 
applied, and with relative accuracy of 0{Aqcd/i^q) in the kinematic region where 
non-local effective field theories can be applied. These expressions agree, within 
this precision, with those obtained from the hadronic result using the layer-function 
approximation plus Euler-McLaurin expansion. Very good numerical agreement for 
the moments is obtained between the exact result and the result using effective field 
theories. 
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A Conyentions and notation 



1 Introduction 



Asymptotic freedom can be seen as the first example of factorization between high and 
low energies, since it dictates that Green functions at high Euclidean energies (Q^) can 
be described by perturbation theory up to corrections suppressed by powers of Aqcd over 
Q. Therefore, the use of the operator product expansion (OPE) in processes where the 
relevant momentum scale is large and Euclidean is safe. This is quite restrictive, since, 
in most of the cases, it can only be tested with experiment through dispersion relations, 
which involve measurements up to arbitrarily high energies. To avoid this problem what 
one usually does is to directly apply the same perturbative factorization techniques to 
observables living in the Minkowski regime. In practice this means to perform the analytic 
continuation of approximate perturbative results obtained in the Euclidean region to the 
Minkowski region. Nevertheless, such calculations do not come from first principles. This 
affects the OPE and effective field theories that are built using perturbative factorization 
techniques aiming to factorize high from low energies. This problem is usually stated as 
duality violations. We will follow here the definition of jl] for duality violations. 

One can quantify the discrepancy between the exact result and the one using per- 
turbative factorization in the large limit of QCD [21 . In this case one finds a clear 
discrepancy between both results in the physical cut of the Green functions, where one has 
infinitely narrow resonances on the one hand and an smooth function on the other. This 
can be further quantified in the 't Hooft model [3], which we will consider in what follows^. 
For this model, when one considers some inclusive quantities like the total heavy meson 
or tau decay rate, the discrepancies between the hadronic and OPE-like result (using the 
OPE for the tau decay and HQET for the inclusive heavy meson decay |4f 5j ) appear to 
be quite suppressed. 

On the other hand one may also study more exclusive quantities like the differential 
cross section of the electron-meson scattering going to electron-|-anything: eM eX 
(deep inelastic scattering), the differential semi-leptonic inclusive decays of heavy mesons: 
Hq — i> Xlu, or e~^e~ — > light hadrons. Indeed one would expect that the magnitude of 
the duality violations for these quantities would be larger, since they are more exclusive 
observables. Some of them have been already studied in the literature |6|7j . like deep 
inelastic scattering or e+e~ — > light hadrons. One finds that the violation of duality is 
maximal but that, if one makes some sort of smearing of the hadronic result, the partonic 
results are recovered at leading order. Nevertheless, for these observables, the partonic 
computation is performed at a diagrammatic level, which makes difficult to go beyond 
the leading order partonic result. In part this is so because one has to deal with jets (very 
energetic final states with relatively small invariant mass) in the final state. Moreover it 
is also difficult to quantify the error made by the smearing procedure, since the smeared 
function does not actually correspond to the differential cross section or decay anymore. 

^One may believe that the large Nc limit sets a kind of upper bound on the duality violations. In the 
real case, the existence of finite decay widths is expected to smooth the duality violations. However, at 
present, it is not possible to quantify this effect. In any case, this does not mean that one can use the 
results of the two-dimensional model as upper bounds to the four-dimensional case with finite Nc. 
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At this respect, there have been recent developments in order to apply effective field 
theories with perturbative factorization to jet physics in four dimensions. This is a rapidly 
evolving field '8'9"i n|lipi2j and the effective field theory has been called soft-collinear 
effective theory (SCET). The use of effective field theories with perturbative factorization 
may allow for a comparison between partonic and hadronic results in a more systematic 
way, beyond the leading order partonic result, and to set up a right framework on which to 
quantify the quark-hadron duality violations. Nevertheless, several questions still remain 
open in SCET, like what the modes of the theory are, or what the optimal formulation of 
this theory could be. The standard formulations of SCET involve the existence of a large 
number of modes, and one can never be sure that the set is complete. For instance, in Ref. 
fr2\ it has been argued that there could be some extra modes called messenger. There is no 
consensus on this issue, though, and in Refs. |13|14j it is claimed that there is no need for 
such modes. Therefore, it is evident that the study of SCET is interesting on its own and 
the application of SCET in a controlled setup may help to better understand the structure 
of the effective theory. Obviously, QCDi+i provides this controlled setup. This will be one 
of the main subjects of this paper. On the side of the optimal formulation of SCET, we 
would like to incorporate the advantages of the light-front quantization frame J3] and the 
associated Hamiltonian-like formulation. This provides non-trivial information, since it 
allows to relate the correlators that appear in the effective theory with the wave-function 
of the bound state. It also avoids to perform complicated diagrammatic computations and 
resummation via Dyson- Schwinger equations to obtain the matrix elements and vertices 
(as it was done in Refs. |6|7j ) . Moreover, working in the light-cone gauge is convenient 
to make the theory effectively abelian in 1-1-1 dimensions. 

The specific observable we will consider in this paper to illustrate the discussion will be 
the differential semileptonic inclusive decays of heavy mesons: Hq Xlv. We will only 
consider the kinematical situation when the invariant mass square of the jet, P|-, is much 
larger than ^^qd ~ 13"^ {P^ is the strong coupling, which has square mass dimensions in 
D = 1-1-1). In this situation we will see that one of the modes of SCET, the hard-coUinear, 
is not a dynamical field and can be integrated out, at least in the light-front frame. The 
final effective theory becomes equal to HQET plus an imaginary vertex. This imaginary 
vertex is local in "time" (in the light-front quantization frame), can be computed order 
by order in perturbation theory, and is able to describe the differential decay rate (more 
precisely, the moments). 

We structure the paper as follows. In sec. 2 we analyze QCDi+i in the light front and 
compute the transition matrix elements. In sec. 3 we compute the hadronic differential 
decay rate and moments. In sec. 4, we work out SCET in two dimensions and compute 
the differential decay rate and moments at tree level. In sec. 5, we develop an alternative 
effective theory without hard-collinear fields and compute the differential decay rate and 
moments at one loop. In the Appendix we set up the notation and conventions. 
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2 QCDi+i in the light front 

In D = 1 + 1, the QCD Lagrangian is given by 

i 

where Z)^ = <9^ + and i labels the flavor. 

One can perform the quantization in a frame different from the equal-time frame. In 
particular it is possible to choose the quantization frame at x'^ = constant, which would 
play the role of time in this case. The role of the energy is played by the conjugated 
variable P~ . The other variables: P+ (and P± in four dimensions) are kinematical. For 
instance, the component of an hadron behaves in "free" -particle way, 

^i^ = E^^ (2) 

i 

where the sum extends over all the partonic components of the bound state. This allows 
to define the variable "x", which measures the fraction of P^ momentum carried by a 
given parton. 

The notation for the components of the gluon field is 

A+ = n+-A, A- = n_-A. (3) 

The usual quantization gauge is = 0. In this situation, the fields ip- = ^-i' (Joy 

the definition of A+ /_ see the appendix) and A~ are non-dynamical and can be integrated 
out from the theory (they are constraints)^. The resulting Lagrangian reads (■?/'+ = A^ip) 

+ / dy-4^t''^,+ {x-,x-^)\x- -y-\ij]^eij,+{y-,x-^) , (4) 

where we have defined 

f -1 , a; < 0, 
e{x) = J , X = 0, (5) 
[ 1 , a; > 0. 

The representation of the quarks in terms of free fields in the light-cone quantization 
frame reads 

Mx) = 1 (a(p)e-^^^ + b\p)e^^^) , (6) 

^One should not forget that there is another constraint, the Gauss law, that restricts the Hilbert space 
of physical states to those which are singlet under gauge transformations. See for instance where 
one can also find a quantization in the path integral formulation. 
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and the anticommuting relations are 

{a(p), a\q)} = {b{p), b^{q)} = 2(27r)5(p+ - g+) , (7) 
{a{p),b\q)} = {b{p),a\q)} = 0. 

Once we have the Lagrangian we can construct the Hamiltonian (in the hght-cone 
frame) 
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dx dy '^pl_^t''^pi+{x ,x+)\x -y l^/j^^fi/jj+iy 



By solving the eigenstate equation (taking into account the constraints and where n 
schematically labels the quantum numbers of the bound state) 

p-\n) = P^\n), (9) 

one obtains the basis of states on which the Hilbert space of physical states can be spanned. 
Here we will focus on the meson sector of the Hilbert space and we will generically 
label the state as \ij;n), where i labels the flavor of the valence quark, j labels the 
flavor of the valence antiquark and n labels the excitation of the bound state. The 
solution to Eq. © can be obtained from the large Nc limit solutions within a systematic 
expansion in 1/Nc using standard time-independent quantum perturbation theory. It has 
the following structure (the momentum of the bound state will not be displayed explicitly 
unless necessary) 

\tj;n) = |u;n)(°) (10) 



+0 



m,n' k 
1 



p{0)- _ p(0)- _ p(0)- 



where the second term in the expression is l/y/N^ suppressed. Here we have used the 
fact that, at order l/^/Wl, P~ only connects neighboring sectors (n-mesons — * n ± 1- 
mesons), becoming an almost diagonal infinite dimensional matrix. \ij]n)'^^^ represents 
the eigenstate solution to Eq. © in the large Nc limit, and Pn^ the associated eigenvalue 
(we do not explicitely display the flavor content of Pi°^ except in cases where it can 
produce confusion). In this limit the sectors with fixed number of quarks and antiquarks 
are conserved and consequently the number of mesons. Therefore, the bound state can 
be represented in the following way 
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where a is the color index, 0^ is the solution to the 't Hooft equation, which will be 
reviewed in the next section, and the state is normalized as 

W (zj; ml^'j'; n) ^ = 2vr2PW+5_5,,5,,,5(pW+ - pW+) . (12) 

The fact that the number of particles is quasi-conserved makes possible to formulate 
the theory along similar lines of how is done in pNRQCD (for a review see pZ|), where the 
wave function (the 't Hooft wave function in our case) is promoted to the status of being 
the field representing the bound state. We will not pursue this line of research further in 
this paper but we expect to come back to this issue in the future. 



2.1 The 't Hooft equation 

By applying the operator P~ to its eigenstate |n) at leading order in one obtains 
the 't Hooft equation 

^ + 7^ C(^) - / dy<^i{y)V- , (13) 

X \-x) Jo \V-^) 

where M„ is the bound state mass, = x = p'^ /P^, with p"*" being the momentum 
of the quark i, and P stands for Cauchy's Principal Part^. The renormalized mass is given 
by m^jij = — P'^. The principal value prescription serves to regulate the integrand 
singularity, which originates in the infrared divergence of the gluon propagator. This 
equation has a discrete spectrum of eigenvalues that increase approximately linearly for 
large n, and the wave functions vanish at the boundaries with the asymptotic behavior 

^^(x)^x^', x^O, (15) 

where 

rn^ R + T^PiCotnPi = , (16) 
and similarly for x — > 1. The 't Hooft wavefunction are chosen to be normalized to unity 

dx(l)^*{x)(l)%{x)=6nm. (17) 

We will have to consider high excitations of mesons when considering the decay rate of 
the heavy meson. In the asymptotic limit n — oo one can obtain analytic expressions 
both for the masses 

M^~n7r2/32, (18) 



■^One can use the following representation of this distribution 

(x - y)2 2 



P-J— = -ll dz|z|e^(--^)^ (14) 
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and for the meson wave functions 

(pliix) ^ V2sm{mTx) . (19) 

Actually, a more detailed study of the 't Hooft equation has been performed in the large n 
limit using semi classical (WKB approximation) techniques in Ref. ^Hl- In this reference 
the layer function was defined (see also [7j): 

0.(0= Iim0jf(e/M„2), (20) 

n—>oo 

for finite ^. This function is the solution of the equation 

m = -f^UO - d^'MOPj^rr^, ' (21) 

and the following equalities can be obtained 

rd^^ = n^, rdmO = P^rn,, (22) 

Jo S "^i Jo 

which we will need in the following sections. 
2.2 Transition matrix elements 

We are now in the position to compute the transition matrix elements due to an arbitrary 
current: 

{cs;m\ipcTQ\Qs;n) , (23) 

where F represents a generic Dirac matrix. We anticipate the notation that we will use 
for the heavy meson decay: Q represents the field of the heavy quark as well as its flavor, 
s the flavor of the spectator quark and c the flavor of the hard-coUinear quark. We will 
restrict ourselves to the kinematical situation relevant for the semileptonic heavy meson 
decay. This means that < and P~ < P^ . 

We only aim to obtain the matrix element ()23|1 at leading order in l/Nc- Nevertheless, 
this does not mean that we can just work with the leading order solution to the bound 
states. As we will see, we will also need the l/y/Wl corrections to the bound state. The 
contribution to the matrix element ()23|) can be split into two contributions. We distinguish 
the contributions to the current according to whether they come from "diagonal" or "off- 
diagonal" terms, which we show in Fig. ^ The diagonal term directly connects the 
current to the leading 0{1/N^) term of the bound state: 

{cs-m\'iPj:Q\Qs-n)Uiag. = (°)(cs; m|4rQ|Qs; . (24) 

This term is of 0{1/N^) and it is produced from terms of the type ipc^Q ~ 6c0.q + ■ ■ • . 
In a way they change the flavor of the bound state from "heavy" to "hard-coUinear". 
Nevertheless, there is another possibility: tpc^Q ~ bcbq + ■ ■ ■ , which can be understood as 
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Figure 1: Contributions to the hadronic matrix elements of the current. The first figure 
corresponds to the "diagonal" contribution to the matrix element, Eq. \24^ . The second 
and third figures correspond to the " off- diagonal" terms, Eq. \25\) . The ® represents the 
current, and the gluon exchange the effective four-fermion interaction in Eq. 
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the creation (annihilation) of a new bound state. This possibihty does not have overlap 
with the leading order term in the expansion of the bound state but it does with the 
1/v^iVc one. Whereas the matrix element connecting the one meson sector with the two 
meson sector is l/^fWl suppressed, the overlap of the two meson state with the current 
is \fWc enhanced. This is why this contribution has to be considered as well at leading 
order in l/N^. We define 

9.(?.7v)P+ p(0)- ^^^^ 



X (0|4rQ|gc; nV°)(°)(gc; n'\^^\cs; m\p-\Qs; n)(°) 



)(0)- 



p(0)- 



A good thing of working this way is that, once {Qc;n'\^'^^ {cs;m\P^\Qs;n)^^^ has 
been computed, it can be used for any current. The total result for the matrix element 
at leading order in 1/Nc then reads 

{cs;m\ijj:Q\Qs;n) = {cs;m\ipcrQ\Qs;n)\diag. + {cs;m\i)j:Q\Qs;n)\off-diag. ■ (26) 

We are now in the position to apply the above discussion to some specific currents. 
We display the results below {x = = (P+ - P+)/P+). 

{cs]m\iicl^Q\Qs;n) = 2{cs]m\ii\^cQ+\Qs;n) 
-1 



(27) 



= 2P+(l-a;) 



dzct>f{x+(l-x)z)^Z{^) 



-x^(3^ 



1 



^0 



dudz- 



'z)G{u-q^ 



^f{x + {l-x)z)-4>^'{xu)) 



x{i-u) + {i-x)zy 



„1 oo 

G{u; q') = dvJ2 

•^0 n'=0 



f - Ml 



{cs]m\ipcn Q\Qs;n) = 2{cs;m\ 



2mQ 



gme /•^^ >f(x+(l- 
L (^+(1- 



id+ 



(28) 
(29) 



x)z)z P+ 



X 



1 rl rl 



^0 ^0 



dydtdz 



Qci 



{til - x) + {\ - z)x)- 



-X<\>f{x-r{\-x)t)-C^f{xz)). 



Jo 



(x + (i-x)^)C(^) 



mQ(l - x) ^ rUc 
X + (1 — x)z z 



(30) 



X 



n' odd 

dydtdz ,5' ^^A^^^y-'l'L (^(^ + (1 - ^)t) - . 



^0 JO 



{t{l - x) + [l - z)xf 
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The "+" component of the vector current was aheady computed in Ref. [?]. The com- 
putation of the rest of the matrix elements had to wait to Ref. P^, but we disagree 
with their results for the "-" component of the vector current and for the scalar current. 
More recently, expressions for the transition matrix elements of the current have also been 
worked out in Ref. [20, • We find that our expressions are more compact than those. In 
any case we have not been able to check the agreement with those. In Ref. j^I] the matrix 
elements have also been considered using similar techniques to ours. Nevertheless, they 
consider different kinematics, which makes difficult the comparison with their results. 

Finally, we would like to stress that current conservation imposes strong constraints 
on the form of the currents. The following equalities have to be fulfilled between the 
different matrix elements 

q^{cs]m\^pcYQ\Qs]n) = {mq - mc){cs;m\i)cQ\Qs;n) , (31) 

where = — P^^. For the case in which the hard-collinear and the heavy quark 
correspond to the same particle, and taking the limit 0, we obtain the equahty 

(first obtained in Ref. |22j ) 

(«m.M^A,«)V)^4^. (32) 

^Ir Jo ^Ir 
2.3 The static limit 

Since in this paper we will study the differential decay rate of a heavy meson, it is 
convenient to consider the specific case on which one of the quarks is very heavy (the 
static limit). If we redefine the heavy quark field, 

Q+ = e-^"^'3^-^Q+„ , (33) 

where in the infinite mass limit one can use {p = rnqv + k) 

^«.(A:)e-'^% (34) 

at leading order in l/mq the Lagrangian reads 

i i 

+^Y. I ^^l/"^l+^>i+(a;",a;+)|a;--rl^,V>,+(2/",a;+) 

dy~4+t''A+{x~,x^)\x- -y-\Q\eQ^Sy-.^'')- (35) 
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In the same way we can obtain the Hamiltonian in the static hmit (at leading order in 
the l/mq expansion) 

Prtatic = -i^^ f dx-dy-i{jl4x-,x^)e{x- -y-)ilJi+{y-,x^)-Ql^X-^d-)Q+^ 
-^Y. [ ^^2/"^l+^>«+(a:-,x+)|x--2/-|^]+t>,+(y-,x+) 

"yE/ dy-4+t''i^i+{^-.^^)\^- -y-\Q\.t''Q+Ay-.^^)- (36) 

Once this Hamiltonian is applied to mesonic states, one obtains the 't Hooft equation in 
the static limit for which one can find a thorough study in Ref. j2Sl- In this limit, one 

works with the function \l/^(t) = ^^0^* ( 1 — — ) , where t = [1 — x)mQ, and considers 

its static limit, which is described by the following equation (e„ = M„ — ttiq, neglecting 
l/mq corrections): 

rri^ - 3'^ ■ t ■ 3'^ f°° (s) 

enKit) = -^^Kit) + -Kit) - ^ "^'JT^^ ■ (2^) 

The quantities that will be needed in the following are expectation values of the variable 
t in the static limit, defined in terms of the heavy meson wave function, 

POO 

(f)^ / dtiKit))"t' ■ (38) 

The numerical computation of these expectation values can be cross-checked with the 
help of static limit sum rules [221 ^^ch as 

(t) = e„ , (39) 



{t') = l{tf-l{ml-P^) . (40) 

The above expressions will be useful in the computation of the OPE expansions of the 
differential decay rate moments, see sec. 13.31 



3 Semileptonic differential decay rate 
3.1 Kinematics 

We consider here the semileptonic heavy meson decay: Hq — > XJJb, where Hq represents 
a bound state made of a heavy quark Q and a light (spectator) quark s (by default we 
will perform the numerical analysis for the ground state but the formulas hold for any 



12 



state). Xc represents any hadronic final state with c (hard-colhnear) flavour content and 
la^h represent massless leptons. We will consider the situation on which the spectator, 
ips, and hard-coUinear, ipc, quarks have different flavour in order to avoid annihilation 
and Pauli interference terms. This decay has already been studied in the past. We will 
follow here the work of Bigi et al. |1] . The authors considered the flavour changing weak 
interaction 

/:l.^ = --^^Pcl,QlaYh. (41) 
The total decay width can be written as 

= ^ 1 ^,9iq^)9{q-)Imn,^{q)ImT'^^{q) , (42) 
where 11^,^ (x) and T^,y(x) are deflned as 

n^,(x) = t (o|T {UxhM^) WhMo)} |o) , (43) 



T^'^ix) = I {Hq\T { Q{x)YM^) ^c(0)7'^Q(0)} \Hq) , (44) 
and their Fourier transform as 

n^.(g) = j d^xe'^'^U^.ix), T^-'iq) = j d'xe-''^^T^''{x) . (45) 
The imaginary part of Il^^{q) reads 

{q~^S{q~) (+,+) component, 
g-5(g+) (-, -) component, (46) 
otherwise . 

We notice that the result Eq. ()46p is the same in the equal-time or in the light-front 
formalism. This is so in the massless case. Once masses are included the situation becomes 
more complicated. One may wonder how it is possible to obtain the term proportional to 
Slq'^) and finite q^ = + = m^/ P^ + ml/ P^ in the leptonic correlator for massless 
leptons. This is somewhat amusing if one works in the light-front quantization frame. 
Naively one would expect that q~ is always zero if the masses of the leptons are zero. 
Nevertheless, one may obtain /mil = q~S{q^) (which is necessary to restore parity 
invariance) by working with flnite masses for the leptons and taking the massless limit at 
the very end. Being more precise, this contribution appears from very high P^^, scaling 

like P+, ~ P^,Jml, with P+, flnite. 

We can see that two terms are generated for the differential decay rate (where Mhq 
is the mass of the Hq meson): 

^-W^""^"""^'''" '''' 
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with X = /Mhq > 0, and 



^xImT++{q-,q+ = 0) (45 



dx 2(47r 

with X = q~ /Mhq > 0. The total decay width then reads 



^H^ = I d4-JZ- + -JZ-] ■ (49) 



Q \ dec dx 
The procedure of Ref. [Ij was to assume that 

dr(+) _ rfr(-) ^ 1 c/r 

dx dx 2 (ix 



(50) 



are equal by parity symmetry and only to compute That both terms are equal is 

indeed highly non-trivial due to the fact that the gauge fixing A'^ = and working in 
the light-front quantization frame breaks the explicit invariance under parity. We will 
here explicitly compute ^^-^ and compare with We then have to compute the 

differential decay rate (in Ref. Q only the total decay rate was considered). As we will 
see, its computation is highly non trivial and requires to take the massless limit for the 
hard-coUinear quark in a careful way. 




Figure 2: Decay of the heavy meson Hq to the meson \cs; n) and the fictitious cj) particle. 



Eq. (j46j) shows that, at the practical level, the interaction could be simulated by 
a massless particle. Therefore, there are some kinematical similarities with a kind of 
h X(Pj decay. The effective interaction would read 0] 

Cak = --^^cl^Qe^^m) + {h.c.) , (51) 
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where represents a fictitious pseudoscalar particle. The process can be seen in Fig. |21 
Since the pseudoscalar particle is real, we can parametrize its two-momentum q as (for 
definiteness we set the kinematics relevant for the computation of for things 

work analogously) 

(f = q = q =0, (52) 

where, by momentum conservation, the momentum of the final hadronic state Px = Pn 
reads 

P°=^(2-a:), Pk = -q\ Pl = Ml^{l-x). (53) 
In light-cone coordinates 

P^=PI-P\ = Mh^, P+ = P0 +Pi = M^Jl-x), (54) 
so we see that in the endpoint region x — > 1, the factor 



is small. It will play the role of one of the SCET expansion parameter, and leads to 
appropriate scalings for the momentum of the final meson state: 

P^~l, P+~P, (56) 

which behaves as a coUinear jet. 

The heavy quark mass will be considered to be a large parameter (equivalent to Q 
in jet physics or deep inelastic scattering). We can distinguish at least three kinematical 
regimes (we use the names SCETI and SCETII for an easier comparison with the notation 
used in effective field theories, see next sections and Ref. 



a) OPE; P| = Ml = M^Jl - x) ~ m| > mgAQco; n ~ ml/g\ 

b) SCETI; P| = Ml = M^Jl - x) ~ mQAQcc; n ~ rriQ/g, 

c) SCETII; P| = Ml = M^Jl - x) ~ Kl^^- n ~ 1. 

One usually refers to situation c) as the most exclusive and a) as the less exclusive 
one. Here we would like to stress that in the large Nc limit (irrespectively of the num- 
ber of spatial-time dimensions) all the three situations are equally exclusive, since they 
correspond to only one physical hadronic final state, n represents the principal quantum 
number of the hadronic excitation (we are having in mind a linear Regge behavior). In- 
deed, the jet multiplicity of the hadronic final state is not well represented in the large 
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We have two independent expansion parameters: A = and A = \/\ — x. We 

note that A << 1 is always fulfilled. With respect A, we will restrict ourselves to the 
situation where we are either in the OPE or SCETI situation. 

Let us note at this stage that we are actually using the opposite kinematic condition 
to the one used in Ref. There is a reason for that. In our kinematics, the "energy" 
Px of the hadronic jet is much larger than Aqcd- Therefore, in the "time"-axis, x"*", the 
interaction takes place at very short times and can be considered local. This is what it 
will allow us to write the interaction as a local term (in "time", i.e. x^) when we try to 
represent the process by means of effective field theories later on. See Fig. E| This will 
also allow us to write the matrix elements in terms of the wave-function of the bound 
state. We will elaborate on this in sees. 4 and 5. 



3.2 Differential decay rate: hadronic computation 

We can use the spectral decomposition to relate ImT with the transition matrix elements 
of the currents we computed in sec. 12.21 We obtain (we have already restricted to g"*", 
g- >0) 

+ a - Pm-<1- + P-H, - Pn) 

X (ifQ|Q(0)7>c(0) |cs; n, P+) (cs; n, P+|V5c(0)7^Q(0) li/g) . (57) 

The expression for the differential decay rate then reads (we work in the rest frame of 
the heavy meson with = P^^ = = P" and x = 1 - P+ZP^ = 1 - Ml/Ml^) 



dx 



327r 



X 



X] 



{n-P^lMOh-QmHg) 



S P, 



p: 



{51 



where the matrix element can be read from Eq. ()29|1 taking the limit ^ 0. We 
notice that the differential decay rate consists of a sum over deltas at the position of the 
resonances and, therefore, cannot be obtained from perturbative-like computations. 

We could also do the computation with the kinematics g+ = 0, g~ = xMhq, along 
the lines of Ref. ^ . This is the " spatial" component of the momentum. In this case we 



would obtain {P^ = Pj 



Hr 



P+ and X = 1 - P-/P^^ 



dx 



r<2 

\j \ — \ 



X 



M„<Mh, 



(n;P+|4(0)7+g(0)|iJQ) 



6 P, 



Hq 



1 - Ml/Ml^ 

Pn - 



(59) 



Q 



Note that in this case we have to compute the matrix elements in the limit g^, P^ 
0, which considerably simplifies the computation and one obtains 



-P^ 



dx 



E 

M„<Aff, 



X 



1 



Ml 



(60) 
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for the differential decay rate, where r„ is 

r, 



dz(l)''^'{z)(l)HQiz) 





(61) 



In principle, this expression should be equal to Eq. ()58|1 for all x. This implies the 
following remarkable identity among matrix elements (x„ = 1 



mgnic ^^(pHQ {Xn + (1 - Xn)z)<p''^{z) 



{P^,)- JO 



n'=0 



1 rl /-l 



^0 ^0 



dydtdz- 



[Xn + (1 - Xn)z)z 



(62) 



{t{l - Xn) + {I - z)Xny 
Xn)t) - (pHgiXuZ)) . 



In the left-hand-side (LHS) of the equality only the "diagonal" term of the "+" current 
contributes. For the right-hand-side (RHS), the first term is the "diagonal" contribution 
of the "-" current and the second term is the "off-diagonal" one. Note that, in principle, 
we cannot fix the relative sign between both matrix elements. Note also that the equality 
Eq. ()62j) provides different information than Eq. (jHT|) . since it relates matrix elements 
with different "x" (P^^ - P+ = in the LHS of the equation and P^^ - Pn = xMhq 
in the RHS of the equation). We have not been able to find a general analytic proof of 
these remarkable identities, though we have been able to do some partial checks, either 
when we have considered moments, or by using the layer functions for the final state (this 
implicitly assumes that we are working with a final state with a large quantum number n) . 
In those cases we have been able to perform a comparison within an expansion in l/mq 
and check the low order terms in this expansion. Irrespectively of the above, we have been 
able to check the equality (jU^ numerically to a level below the 1 % using the numerical 
solution to the 't Hooft equation obtained from the Brower-Spence-Weis improvement of 
the Multhopp technique ^Hj. We have used two set of values for the masses of the quarks: 
TTiQ = 15/9, rric = 10/9, = 0.56/3 and rriQ = 10/3, rric = /?, = /?. We show the 
comparison in Tables Q and |21 This agreement is quite remarkable if we take into account 
that the support functions in both integrals are quite different (see Fig. El), specially for 
the second set of parameters. Note that we do not consider the second term of the RHS 
of Eq. ()62j) in the plot. This term appears to be a correction compared with the first 
one and vanishes in the limit mq oo. This is illustrated in Tables Q and |21 where this 
second term appears to be smaller for larger values of the heavy quark mass. This points 
to the fact that their scaling may go like ~ (3'^/mQ and that there are no terms of the 
type ~ ml/rriQ. The convergence of this second term is very slow. One has to sum over 
a very large number of states to converge to the final value. We illustrate this problem in 
Fig. EJ where the sum is over 100 states. 
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n 


LHS 


RHS ("diag" term) 


Rel. Err. 


RHS 


Rel. Err. 





0.96433 


0.97144 


7 10-3 


0.96627 


2 10-3 


1 


0.25999 


0.26175 


6 10-3 


0.26025 


1 10-3 


2 


0.04432 


0.04473 


9 10-3 


0.04445 


3 10-3 


3 


0.02389 


0.02405 


7 10-3 


0.02391 


1 10-3 


4 


0.00543 


0.00538 


7 10-3 


0.00540 


3 10-3 



Table 1: Values for the matrix elements as defined in Eq. h6^] . The first column 
corresponds to the principal quantum number. The second column corresponds to the left- 
hand side of the equality. The third column corresponds to the first term of right-hand 
side of the equality (the "diagonal" term). The fourth column to the relative difference 
between the second and third column. The fifth column corresponds to the right-hand side 
of Eq. / I6'ij)) and the last column to the relative difference between the left and right-hand 
side of Eq. h6^) . In order to ease the comparison with the results of Lebed and Uraltsev 
IB/, we take mq = 15/3, rric = 10(3 and nig = 0.56/3. 



n 


LHS 


RHS ("diag" term) 


Rel. Err. 


RHS 


Rel. Err. 





0.46946 


0.42992 


9 10-2 


0.46493 


8 10-3 


1 


0.61406 


0.62957 


2 10-2 


0.61537 


2 10-3 


2 


0.49594 


0.51617 


4 10-2 


0.49821 


3 10-3 


3 


0.32820 


0.34773 


6 10-2 


0.32985 


5 10-3 


4 


0.18571 


0.19712 


6 10-2 


0.18694 


6 10-3 



Table 2: As in TableUl with the values mq = 10/3, nic = [3, rus = p. 
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Figure 3: Plot of the integrands of Eq. h6^] . We have taken the values n = 3 and 
rric = 10/3, rris = 0.56/5 and uiq = 15/5 in the first figure and = 1/5, = 1/5 and 
mq = 10/5 in the second figure. The dashed red line corresponds to the integrand of the 
left-hand side of the equality. The solid blue line corresponds to the integrand of the first 
term (the second term is subleading in l/rnq and it is not considered in this plot) of the 
right-hand side of the equality. In the first figure the solid blue line diverges (although in 
an integrable manner) for z —>■ but it cannot be seen with the resolution of the plot. 
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Figure 4: Analysis of the (off-diagonal) second term in the right-hand side of Eq. / I6'ij)) . 
We have taken the values n = 3 and = (3, = (3 and tuq = 10(3. The dashed blue 
line corresponds to the difference between the left and right hand side of the equality if 
the off-diagonal term is neglected. The solid red line represents the contribution of the 
off- diagonal term as a function of the number of intermediate states added to the sum. 
For n' larger than 80 the numerical stability of the computation is doubtful. 



We are also able to compare with the numerical evaluation of the matrix elements 
(the LHS of Eq. (jHSI)) performed in Ref. jSj. We have checked that our results agree with 
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theirs within the expected numerical uncertainties. 

We would also like to remark that the RHS of Eq. (jU^ can be understood as a function 
of X, which for equal to the LHS. Therefore, it provides with a definition of 

a continuous function in x. This will be relevant later on when trying to connect with 
computations using effective field theories. 

So far these expressions are exact. At this stage we can perform an expansion in l/mg 
and consider the large n limit (therefore our result will hold for the OPE or SCETI region 
but not for the SCETII kinematical situation). At lowest order in those expansions, and 
using the properties of the layer function defined in sec. 12.11 we obtain for the "diagonal" 
term of the matrix element (we also include the subleading corrections in mc/mQ, which 
can also be reliably computed with the layer function) 

m,m. rju,{x + {l-x)z)4>-{z)^rn^ T ci.^^C f^V) 



Jo {x + {l-x)z)z Mi^Jx y y-x^'^yi-x 

mgm, ( <Phq{x) di 1 l-x 

L y0c(o+0H,(a:)^^y^ dim) 



x^ 



Ml JO 



d^MO ]+■■ 



rriQ 0Hg(x) / ml 1 ml 

vrp 1 + —5 — ■ — —5 h 



Ml^ X ^\ Ml^cl>H^{x) xMl 
by defining 

y-x ^ J_ 
1 - X ~ M2 



(64) 



and expanding in 1 — x. 

We would like to stress that both the LHS and RHS of Eq. (jU^ can be defined for any 
X and not only for x = Xn- Nevertheless, they correspond to the physical matrix element 
only for x = x„. On the other hand the RHS of Eq. ()63|) provides with an interpolating 
function for the matrix elements at different x„, which is independent of the dynamics of 
the final state 0„. We plot this function in Fig. El We can see that the next-to- leading 
order (NLO) is a correction compared with the leading order (LO), for all values of x 
and the hard-collinear mass we consider. For instance, if we take x = 0.5, we obtain 
0.050022 = 0.0470693 + 0.0029527, where the first term is the LO resuh and the second 
one the NLO correction. It should be noticed that most of the contribution to the NLO 
result comes from the derivative term of the wave-function. We will discuss further this 
issue in sec. 13.31 when we consider the moments. 

So far we have only considered the "diagonal" term of the hadronic matrix element. 
For the "off-diagonal" term, we are, at present, not able to give approximated analytic 
results, even in the large n and large mg limit. Nevertheless, we have some hints about 
its analytic form in those limits. They come from two sources: a) the expressions for 
the moments from the hadronic computation which are already available, for = 0, 1 
and 2, within an expansion in l/mg and b) the results from the effective theory with 
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Figure 5: The solid red line is the LO expression of the RHS of Eq. h6^) . The dashed 
blue line is the NLO result. We take mq = 10/5 and nig = [3, whereas rric = P, 2/5 and 3/3 
for the first, second, and third figures, respectively. 



one-loop accuracy. As we will see, to get agreement with these results, the leading term 
in the large n and large ttlq limit of the "off- diagonal" correction should renormalize the 
masses of the hard-coUinear and heavy quark of the "diagonal" term. In practice, one 
should have (although the absolute sign cannot be obtained, the relative sign with respect 



21 



the "diagonal" term is fixed): 



{t{l - x) + {1 - z)xy 



\-^Hq) „,^o 



Ml, 



X 



/?2 1^^ ^ 



(65) 



Although we were not able to check this equation on an analytic basis, we have been able 
to check it on a numerical basis. We show this comparison in Fig. IHl We can see that the 
LHS and RHS of Eq. converge to the same value as expected. We have also checked 
that if we vary the masses^, the same pattern survives (as far as the heavy quark mass is 
large enough). Note that this term is a correction in the l/mg expansion. 



. 04 



0.03 



. 02 
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Figure 6: Plot of the absolute value of the LHS (dashed line) and RHS (solid line) of 
Eq. i65]) for x = a;„ = 1 — M^/M^^. To improve the numerical accuracy, we use the 
difference between the LHS and the first term of the RHS of Eq. \6^) for the numerical 
value of the LHS of Eq. \65]) . We take the values mq = 10(3, m,, = (3 and nig = (3. 



Overall, we find that the total=" diagonal" +" off-diagonal" matrix element can be writ- 
ten in the following way (up to a global sign) for large n and mq: 



I 




Mfj^x 

Note that the RHS of Eq. (jfifjj) can be understood as a function of x. Strictly speaking 
this expression is singular for x ^ 0. Nevertheless, this effect only shows up for very 

"^One should note though that the computations of the 't Hooft wave function with tachyonic masses 
are problematic at the numerical level. This problem affects the accuracy of the numerical results and is 
more acute if one consider the derivate of the wave- function on a point-to-point basis. 
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small values of x, which are not included in the plots. We can check how well this curve 
compares with the exact hadronic computation for different values of the hard-coUinear 
and heavy quark masses. We show the results in Figs. [3 and |H1 On the one hand we 
plot the hadronic matrix elements: The LHS of Eq. (jH^ and the diagonal term of the 
RHS of Eq. (jn21)- This allows to visualize the difference of working with renormalized 
and not-renormalized masses (actually the only place where this difference is visible is in 
the hard-coUinear mass multiplying the derivative of the 't Hooft wave function of the 
Hq meson). On the other hand we plot the function obtained from the boundary-layer 
approximation, Eq. (jHEI), for all values of x, also with renormalized and not-renormalized 
masses. We can see that the agreement between the layer-function and the hadronic 
result is very good up to very low values of n or, in other words, up to quite near the 
X — 1 limit^. Our results are also quite good up to relatively low values of the heavy 
quark mass. We can also see that the dependence on the hard-collinear mass is very well 
understood with our analytic formula. Moreover, we can also see how the effect of the 
non-diagonal term of the RHS of Eq. (j62p is equivalent to renormalizing the heavy quark 
and hard-collinear mass in the region where we can trust our results^. Overall, we get a 
very consistent picture. 

We find the equality Eq. (jU^ quite remarkable. It implies that the partial decay width 
r„ becomes independent of the final state wave function properties for higher excitations. 
The dependence on the final state only appears through x„ = 1 — M^/Mjj : 



Q,R 



n'p' 1 



47r 



1 + 2 




The differential decay rate then reads 



dx 



X 



1 + 2 





6 



X 



. Xr, 



(68) 



(69) 



M2 

n 

Ml 



^Actually, the agreement is even too good for the x — > 1 limit, where the layer-function approximation, 
in principle, does not apply. One could not ruled out this to be a numerical accident. For instance, if we 



use 



dz(t)n{z)(t>HQ{z) ^ 71-/3 



1 



Q,R, 



(67) 



instead of Eq. H66|l . which is also correct to the order or interest, the agreement is less good in the x ^ 1 
region. On the other hand, Eq. (|67|l incorporates subleading partial effects which may jeopardize the 
agreement. 

^In the effective field theory, this renormalization would be produced by one-loop 0^) corrections. 
We can see that their effects are very tiny for basically all values of x. In the hadronic computation it 
reflects in the fact that the "off-diagonal" effects are also very small. Actually, the basic effect that it is 
seen is the renormalization of the hard-collinear mass that appears in the hard-collinear propagator in 
the effective theory. This is what is to be expected in the SCETI region, in the OPE region the I/tiq 
corrections are too small to be seen by the eye (unless the hard-collinear mass is large enough). 



23 



0.3 



. 1 



0.2 - 



. 4 



0.7 



0.5 



. 6 




X 



0.2 



. 4 



. 6 



. 8 



1 



Figure 7: The solid blue lines represent Eq. / I6'6|) for values of = j3, 2j3 and 3(3. The 
dashed red lines represent Eq. i66]} for values of m^. = (3, 2(3 and 3(3 replacing all the 
renormalized masses by their bare values: rriiji rrii. The dots or squares represent the 
LHS of Eq. for different values ofn. The triangles represent the diagonal term of the 
RHS of Eq. KU^) . In both cases the values of nic are (3, 2(3 and 3(3. We take mg = 10/3 
and rus = (3 in all cases. 



This expression will be suitable to a smoother connection with the computation using 
effective field theories. 

Eq. (jU^ applies to the kinematical situation 1 — x ^ P'^/mg. It includes the leading 
term in an expansion in l/niq and 1/n. Some kinematical l/mq corrections are auto- 
matically included by working with the exact ((>Hq wave function instead of working with 
the strict static limit. The corrections of order ml/mQ, (3'^/mQ have also been included. 
In principle, this expression could be systematically improved by considering corrections 
in l/rrtQ and 1/n. Nevertheless, this would require to know the corrections in 1/n to the 
integrals that appear in our expressions, which at present are not known. We expect to 
study further this issue in the future. 

Finally, we would like to stress that the limit — > has to be taken with care, as it 
is evident from Eq. (j62j) . A naive limit nic — > may lead to wrong results. 



The differential decay rate is not a very well defined object in the large N^.., since it becomes 
either infinity or zero. In particular its comparison with the expressions obtained from 
effective theories that use perturbative factorization is not possible, as we will see in the 



3.3 Moments 
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Figure 8: The solid blue lines represent Eq. / I6'6|) for values of niQ = 10(3, 6(3 and 3(3. 
The dashed red lines represent Eq. KU^) for values of nig = 10(3, 6(3 and 3(3 replacing all 
the renormalized masses by their bare values: rrii^R — ^ mj. The squares or dots represent 
The LHS of Eq. / I6'ij)) for different values of n. The triangles represent the diagonal term 
of the RHS of Eq. KU^) . We take rric = (3 and nis = (3 in all cases. 



next sections: on the one hand one obtains a series of delta terms, whereas on the other 
one gets an smooth function of x''. At this respect one may think that it is better to work 
with moments'^: 



Exact expressions for the moments in terms of hadronic matrix elements can be ob- 

^Real experimental data on semileptonic B meson decays is usually available in terms of moments, 
and therefore so are the corresponding theoretical predictions. In Ref. the differential decay rate 
itself was reconstructed from available experimental information on its moments, allowing thus for a more 
general comparison between theory and experiment. 

®This actually does not cause the problem of quark-hadron duality to vanish, though, as it has been 
emphasized in Ref. for the inclusive decay width. 




(70) 
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tained by using the expressions for obtained in the previous subsection: 



M, 



N 



in 



E 



X. 



N 



M„<Mh. 



Q 



y-^HQ) „,^o 



(X„ + (1 - Xn)z)z 



(71) 



HqJ jo 
1 /.I /.I 



- x„) + (1 - z)x„)^ 



X {(j)HQ {^n + (1 - a;„)t) - (pHn {XnZ)) 



TV 



E 

Mn<MHr 



dy<Pniy)<pHQiy) 



(72) 



where = 1 — M^/Mfj . As we have mentioned previously both expressions for the 
moments yield the same result. 

In order to perform some analytical analysis, it is necessary to be able to compute 
the different matrix elements or, at least, the sum of matrix elements that contributes to 
the moments. In general this is not possible. Nevertheless, for some specific cases, it is 
possible to obtain approximated expressions. If we take Mi from Eq. (f7^ . it corresponds 
to the total decay width. In this case it is possible ^ to obtain a closed analytic expression 
up to 0(1/M^^) suppressed corrections in terms of expectation values of matrix elements 
of the i^Q-meson wave function by using sum rules: 



Ml 



{ml 



Att M, 



Hq 



^ dr 





Q 



M, 



Hq Jo 



E 

Mn>MH, 
1 



Q 



Ho = 77 / (PHr.{^)- 7 , -Ln (73) 

(74) 

4:71 mq 

is the free heavy quark decay rate, and r„ has been defined in Eq. (note that in Eq. 
(I73|) they represent partial decay widths that are not allowed by momentum conservation). 

For 7^ 0, in general, it is not possible to follow the same procedure, since the sum 
rules become divergent. Only for iV = 0, 2, it is also possible to obtain a finite result: 



where 



G"^ mi. - mi 



Mn 




(76) 
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Mo 



G^Mh^ {ml 



m 



2\2 i-l 



An 



dy ,2 ( N 1 



J2 {Ml^-Ml)V^ 



1717 



(77) 



The = case is basically due to probability conservation. It can be noticed that the 
above moments can be written in a more compact way in the following form: 



M. 



N 



4:71 




for AT = 0,1, 2. (75 



For larger than two the integral becomes divergent. 

The above expressions contain some implicit dependence on the heavy quark mass, 
since, so far, we have used the exact i^g-meson wave function. If we perform an explicit 
expansion in l/mq, one obtains, up to ©(l/mg). 



Mn 



Ml 



An 



An 



(t)2-(t2)+/?2 



2m? 



M, 



An 



2mf 



(t) , 3(t)2-3(t2)+3/?2 



1- — + 

mq 



2m2 




(79) 



^0) 



(81) 



where the static limit expectation values are defined in Eq. (|38p . We relegate the numeri- 
cal comparison of these expressions with the exact ones to sec. 15. H since these expressions 
will also be obtained from the effective theory computation. 

The above expressions for the moments have been obtained for low A^. Therefore, 
they correspond, somewhat, to the kinematical regime where the OPE is valid. We may 
consider to use the approximated expression obtained using the properties of the final 
state wave function for large n, i.e. the layer function, for in Eq. (j69p . Therefore, 
we expect the expressions that we will obtain to be also valid for larger values of A^: 
N < MqI (3 (1 — X ~ (3/mQ), up to corrections of order NP'^/mq ~ f3'^/M^. This means 
the kinematical regime where the OPE and SCETI are valid. Note, however, that we 
integrate for all x in the moments. Therefore, this includes contributions from x ~ 1, 
equivalent to final states with n ~ 1 for which the layer-function approximation is not 
valid. A very rough estimate sets the contribution of these states to the moments (for 
A^ ~ 1) of Oll/rriq) or smaller. In any case the fact that we have problems to obtain 
approximate analytic expressions for Eq. ()65|) sets the accuracy of the calculation. The 
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expressions that we obtain for the moments read 



M, 



N 



Att Mfj M 



1 + 2- 
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c,R 



47r M 




X 



(82) 



6 X 



n 



+ 



N 



n- Hq 



2 N 2 

■c,R 



, Ml 
a; - 1 + 



ML 



+ 



where in the second equahty we have reshuffled the NLO correction to the layer function 
in a way that is correct at the accuracy of the calculation and that it will ease some 
intermediate analytic computations, making them more compact. 

By working with moments, which imply an integral over all x, it becomes possible to 
perform a Euler-McLaurin expansion for the sum over n, which, at lowest order, it is just 
equivalent to replace the sum by an integral: — >■ J dn. This replacement allows us 
to make quantitative the comparison between the perturbative and hadronic result. Note 
however that the Euler-McLaurin expansion is an asymptotic expansion. Therefore, it is 
a difficult question to assign an error. Here we will not dwell further on the error used 
by replacing the sum by the integral. To go beyond this approximation would require a 
better knowledge of the properties of the layer function and a systematic procedure to get 
corrections from it, which is relegated for future work. In any case, the result we obtain 
after the smearing reads 



M 



N 



4tt M 



He 



N 




N 



m 



c,R 
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-X 



xm 



Q,R 



x^ 



X + 
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47r M; 



m 



Q,R 



47r M 
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Q,R 



N 



^ dx 



(83) 



^"'^i, (^) 



where in the last equality we have used 



c,R 



Q,R 



which is correct with the accuracy 



of our calculation. Eq. (jHoj) is correct at leading order in the OPE and SCETI kinematic 
region. In the OPE region is correct up to, and including, 0(1/771^) corrections in the 
situation when it is possible to compare with the already known hadronic results [N = 
0, 1, 2). Note that in order to get this agreement it is crucial to " renormalize" the masses 
of the hard-collinear and heavy quark. This renormalization effect can be traced back 
to the "off-diagonal" contribution to the "-" current. We can also give expressions for a 
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general within an expansion in I/itlq. At ©(l/mg) we find the following expression 



9 



^ (2(iV - 2) + 3) {tf + ((AT - 2)(Ar - 3) - 3) {t") 




4^ ^//g ^0 



The expression for M^r, Eq. (jHH|) . also applies to the SCETI region. As we have already 
mentioned, the above expression is correct at leading order in the P'^N/rriQ expansion. 
Our expression also includes the subleading corrections of 0{f3/mQ). Formally, in this 
kinematical regime we could approximate M^r by the following expression (note that 
(pQ (x) should also be expanded in l/mq) 

(^'2 1\/T rn"^ /"l 

M^^-^ - dxe'-^'-'^%ix) (85) 

x|l + 2(l-a:)-Ar^i^-ArI^ + 

where the terms neglected are of relative order which, in principle, we cannot 

claim to have all of them. The reason is that we have not considered terms of the type 
f3'^/mQd'^(j)'jj^{x)/d'^x, which would contribute to the moments at NNLO in the SCETI 
region. Note again the necessity to renormalize the hard-collinear mass. This effect can 
be traced back to the derivative term in the boundary layer approximation of the "off- 
diagonal" term and can be unambiguously identified numerically. The reason this term 
is enhanced is because one has contributions like 

' t/xa;^-2^^k^i!l = _(iv - 2) rfxx^-30^^(x) . (86) 

Finally, we would also like to consider another observable that is usually used in the 
study of the differential heavy meson decay rate. Since, in real life, many times one 
cannot measure over all the spectrum of final particles, one has to introduce a cutoff to 
the inclusive measurement. Therefore, the following observable is usually considered: 

^H^{y) = ^ dx—{x), 0<y<l. (87) 

i Hq Jl-y U-X 

We relegate a numerical analysis of this observable as well as of the moments to sec. 15.11 

^Quite remarkable, we obtain the same expression if we extrapolate Eq. (|78|l to values of N different 
of 0, 1, 2, for which it was originally obtained, and we expand in (1 — x) before doing the integral. 
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SCETI fields 



Hard-coUinear light quark 
Hard-collinear gluon 
Soft light quark 
Soft gluon 



Table 3: Relevant fields in SCETI. 
SCETI modes 



Mode 


p- 


P± 


p+ 


p2 


Hard 


1 


1 


1 


1 


Hard-Collinear 


1 


A 


X' 




Soft 


A2 


A2 


A2 


X' 



^hc 
^hc 

Qs 



Table 4: Relevant momentum configurations (modes) in SCETI. 

4 SCETI: Multimode approach 

We want to describe now the hadronic results obtained in the previous section using 
effective field theories. The aim is to describe the decay of the heavy meson to a bunch (one 
in the large N^) of hadronic particles with invariant moment P|- = Mjj^{l — x) ^ ^qcd- 
The usual procedure, this scale being much larger than Aqcd, is to use perturbative 
computations. Actually, here lies the heart of the problem of quark-hadron duality, since 
we are working in the Minkowskian region and, therefore, near the mass-shell region. 
Nevertheless, we expect that by working with the 't Hooft model we may better visualize 
the problem. 

In this section we first approach the problem adapting the present formulations of 
effective theories for very energetic collinear particles |9|1U|1HT^ . in particular of Ref. 
fll|, to the two dimensional case, and we relegate an alternative approach to the next 
section. In those references, one attempts to explicitly obtain all the modes that one has 
in the theory from perturbation theory. This heavily relies in the concept of threshold 
expansion of Feynman diagrams |24j . 

In this paper we will not exhaustively explore the different modes that may appear in 
the 't Hooft model. We will see that at the order we will work here it will be enough to 
work with hard-collinear and soft modes^", for which we set the notation in Table 121 Our 
aim is to try to see explicitly at which point in this effective field theory derivation one 
approximates the hadronic result by a partonic one. 

First of all let us set the terminology of the different modes. We have already men- 
tioned that the small expansion parameters are A and A. We will formally work in the 

^^Nevertheless, we believe the 't Hooft model provides a nice framework on which to explore which 
modes really appear in SCETI. We expect to pursue this line of research in the future. 
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situation A ~ A ^ 1, which corresponds to the SCETI kinematic region. Nevertheless, 
our resuhs will also be valid in the OPE region where A ~ 1. In terms of A, the mo- 
menta of the different modes scale as shown in Table lU The only difference in D = 1 + 1 
is that there is no perpendicular component p± of the momentum. The momentum is 
decomposed as 

= n+p-^ +p±+ n-p-^ , (88) 

so that a collinear particle is defined as a particle with large light-cone momentum in the 
direction. The fields that are relevant for the effective theory are shown in Table 01 
The field of each particle is decomposed into all the possible modes: ip = ^hc+Vhc+(ls + - ■ ■ ■ 
In practice only a few modes contribute to a given field. For instance, the quark s can 
be approximated by its soft mode: ips — ?s,s- We also use the notation p+ = n+p, and 
the same for other vector components. A field with momentum p varies in position space 
according to the uncertainty principle, 

a:^ = (x+,x^,x-)~( — ,— ,^) . (89) 

p p± p'*~ 

The scaling of the quark and gluon fields in D spacetime dimensions can be first naively 
estimated from the free quark and gluon propagators in position space quantized in the 
equal-time frame^^, 

(o|r {^{x)m |o> = / (01^^-"-^^"^^ ' (90) 



{0\Tmx)A'^{y))\0) 



d'^p 



-g^-" + (1 - a) 



where the gluon propagator is obtained in a general covariant gauge. Using these propa- 
gators, the scaling of the different fields of the effective theory is the following: 

• Soft light quark ~ A-^^^ 

• Soft gluon As ~ A^-2 

• Hard-collinear gluon ~ XiD~A)/2^ ^± _ ^(^-2)72^ _ 

• Hard-collinear quark E,hc ~ A^^"^)/^, rjhc ~ A^/^ 

To obtain the scaling of the hard-collinear quark fields, we have decomposed these fields 
using projection operators 

i^hcix) = ihc{x) + rjhc{x), ^hc{x) = A.tphc{x), r]hc{x) = A+iphcix) . (92) 



^^In the light-cone quantization frame things are more comphcated, in particular for the fermions, 
which have different free propagators. 
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Scaling of the SCETI fields 



Field D = 4 ~ A°] D = 2 ~ A 



^hc 


A 


AO 


Vhc 


A^ 


A 




A° 


A 




A 






A^ 


A3 




A3 


X' 


gA': 


A^ 


A2 



Table 5: Scaling of the effective theory fields in D = A and D = 2 taking into account the 
scaling of the coupling. 



Scaling of the integration element 



Fields 




X 


he 


A- 


-2 


s 


A- 


-4 


hc+s 


A- 


-2 



Table 6: Scaling of the integration element in the effective action. 



One can check that for = 4 one recovers the usual scalings of the effective theory. The 
scalings of the effective theory fields in D = 4 and D = 2 can be seen in Table El One 
can observe that in D = 2 the hard-coUinear gluons have negative scalings. The solution 
to this problem comes by realizing that in D = 2 the strong coupling g has dimensions of 
mass and, actually, sets the scale of Aqcd so 

gr^X\ (93) 

and we can see that the scaling of the hard coUinear gluons times g is positive. 

It can be observed that the hard-coUinear gluon field times g does not scale as the 
corresponding hard-collinear momentum but rather, it is suppressed by powers of A. This 
will lead to the result that hard-collinear interactions are suppressed with respect to the 
corresponding kinetic terms. 

A crucial point to fix the scaling of each term in the effective action is to know the 
scaling of the integration element d'^x, which depends on the fields present on each term. 
The different scalings can be seen in Table 

We are now in the position to write the effective Lagrangian, which we do in the next 
section. 
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4.1 Lagrangian and heavy-to-light current 



We first want to translate to two dimensions the standard procedure to obtain the effective 
Lagrangian. This first means to integrate out the rjhc component of the hard-collinear 
field^^. We will only consider the leading order Lagrangian here. 



\ \ \ 




/ / / 



Figure 9: Allowed vertices by momentum conservation in a theory with only hard-collinear 
and soft modes. Dashed lines can be either quarks or gluons 




Figure 10: Allowed scattering processes. 



The allowed vertices by momentum conservation are drawn in Fig. |H1 Since the gluons 
do not appear as physical particles, we are always faced with diagrams of the sort of those 
shown in Fig. ^1 By power counting we can easily see that the self interactions between 
hard-collinears (in particular those including hard-collinear gluons) are suppressed by 
powers of A. Therefore, at leading order, we only have to consider soft gluons and quarks 
and hard-collinear quarks. The leading order Lagrangian then reads 



£(0) ^ _ ^tr [F^^^sF^"''] + qs,silfsqs,s + Chqet , (95) 

^^This could be considered somewhat strange. If we were working in the hght-front Hamihonian frame 
with hght-cone gauge A'^ = 0, rjhc would correspond to the physical componentent, 'ijjhc,+ , of the field. 
Therefore, we would be integrating out the physical component of the field and keeping the constraint. 
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Figure 11: Tree-level matching of the heavy quark current. 



he - ^tr [{D^A,^ 



c,hc 



2 Sc,hc 



(96) 



In this expression we have kept kinematical subleading corrections proportional to the 
mass of the hard-colhnear. 

We do not exphcitly write the HQET sector of the theory. In practice it will turn out 
more convenient to implicitly keep the heavy quark mass dependence and to expand at 
the end of the calculation. 

The next step in order to apply SCETI to the semileptonic decay is to write the heavy- 
to-light current ipJ^Q (where F = 7^ in our case) in terms of the effective theory fields. 
The emission of a hard-coUinear quark (see Fig. ITTj) by the near on-shell heavy quark 
puts it off shell, and it stays off-shell when subsequent hard-collinear and soft gluons are 
emitted. Therefore the effective current must reproduce these diagrams that are absent 
in the effective theory, and it can be shown that 



QCD 



'i^cT 1 



1 



ilj) - mQ(l - f) 

We must expand the above expression in powers of A. Starting from 

1 



g4hc Qv ■ 



Q= l-T 



il^ - mQ(l - f) 



g4-hc Qv 



and defining 



^0 = 



1 



1 



"^q(1 -|^) 



184 



(97) 



(9^ 



(99) 



where the scaling is like this because the momentum in the covariant derivative is hard- 
collinear, since it is the momentum that flows through the heavy quark line once the 
hard-collinear gluon has been emitted. Expanding in A, one arrives to the following result 



Q 



1 - S^gA^^J^ + SogAhcJ^S.gAncJ^ + 0{\') ) (1 + 0{\)) Q„ . (100) 



In the last term the 0{\) indicates possible contributions coming from the integration of 
the short component of the heavy quark field. 
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We would like here to stress the difference with the situation in 4 dimensions, where 
a hard-coUinear Wilson line appears multiplying the effective current. Actually here this 
also happens but the exponent is suppressed by powers of A (in four dimensions it is only 
suppressed by powers of as). 

This is the first part of the construction of the effective current. The second part 
consists on the matching of the light quark field ipc in terms of the effective theory fields. 
Since the effective current must be constructed such that it reproduces the on-shell matrix 
elements with a current insertion of full QCD, we must add the following interaction term 

£, = V^,rQ5 ^ (el,,, + vUc + <s)j, J = e^-^^^^/rSi? , (101) 

where B = —G/\/2la'^plb in our case, to the Lagrangian Eq. (j94|) . and perform again 
the relevant manipulations (integrating out the small component field rjhc and multipole 
expand), taking now into account the presence of the source term. Now the equations of 
motion lead to ^ 

VcM = -j^ [gAhc+qc,s - j] , (102) 

%L) 

and inserting this in the effective Lagrangian results in the following modified source term 

A- = {iUc + ?l + is9Ahc^^^ J , (103) 

where the derivative in the last term is hard-collinear and acts to the left. Therefore the 
full QCD light quark is matched in the effective theory into 

^ Cm + lis + ql,s9Ahc+-^ . (104) 

Putting everything together results in the following current 

Jqcd{x) = e-*'"«^-' [Oo + Oi + . . .] , (105) 

where the terms are labeled with their relevant order with respect to the dominant term, 
which reads 

Oo = el,,.7°rg. . (106) 

Note that we have neglected the term ql^TQ^, since by our kinematical assumptions, there 
is large momentum transfer to the final state, so the operators of the effective current must 
contain at least one hard-collinear field to contribute to such final states. The scaling of 
each term is Oq ^ X and Oi ~ A^. The final conclusion is that up to gauge invariance 
subtleties, the heavy-to-light QCD current is simply matched to 

jQCDix) ^ e-^"^«''^eI,/.c7°rQ, = e-*'"«^^'Oo , (107) 

at tree level. This is the result that we need to study factorization in heavy-to-light decays 
at leading order. Hard fluctuations would be included in the Wilson coefficients Ck of the 
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operators in the effective current and can be determined by matching calculations so that 
the general structure of the current reads 



Jqcd{x) = e-^'"<3^"5^Cfc(mQ)Ofc . (108) 

k 

Using field redefinitions, one can see that in SCETI in D=4 at leading power the hard- 
collinear and soft degrees of freedom decouple [TU]. In 0=1+1 the same appears to be 
true. The starting point is Eq. (j9l|) . One can see now that redefining the hard-coUinear 
field ^c,hc using a soft Wilson line, 

icMi^) ^ Y,{x)ifhx) , (109) 



dsAs-{x^ + sn_) j , (110) 
and redefining the hard-coUinear gluons in the following way: 

^LW=n(x)4cM(:^), (111) 

leads to factorization of hard-coUinear and soft modes at leading order in the effective 
Lagrangian, 

= + a (^^- - ^) ^ + [(^"4?" " ^"^^ )^] • (112) 

The factorization at the level of hard-coUinear fields is somewhat academic at this stage, 
since they are not going to appear at leading order in the semileptonic decay of the heavy 
meson. Now all the soft-hard coUinear dynamics are encoded in the effective current, that 
at leading order becomes 

= nVrg. . (113) 

In the following subsection we analyze the semileptonic decay using SCETI. 



4.2 Semileptonic differential decay rate 

In this section we show how factorization can be implemented in this process using SCETI. 
In sec. El we have already written the differential decay rate in terms of the imaginary 
part of T . We can now write this hadronic correlator in terms of the SCETI fields. 
The first step consists on the factorization of the hard modes. We have seen that the 
QCD current can be expanded in a series of operators in the effective theory. Then the 
hadronic tensor can be written as 

T-^i I d?xe^^-^^-^^> J2 H,{mQ){HQ\T{Oi{x)0,>m\HQ)^Y.Hk{rnQ)T~ 

k=k'+k" k 

(114) 
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where the effective hadronic tensor at leading order is given by 

T--^ = z I d\e^^-'^+^^^> (ifQ|T{(e,,,,7~Q.)^ (x) {icMl-Qv) (O)}!^^^) • (115) 
After performing the field redefinitions shown in Eq. ()109p the effective tensor reads 

T^^ff = ^ I dV^-^^-'^^^- {Hq\T{QIJx)Y{x)Y\0)Q,AO)}\Hq) (116) 

X (o|r{(7-eSj.(x)(gt7-)/.(o)}|o; 



where we have used the factorization that the (redefined) soft and hard-collinear modes 
hold at leading order at the Lagrangian level. Therefore the correlator can be understood 
as the convolution of the soft and jet function. The heavy quark correlator is explicitely 
soft gauge invariant due to the Y string. So far the computation has been pretty much 
similar to the one in four dimensions (in the four-dimensional case there are also hard- 
collinear strings that here have already been approximated to 1). 
Let us now discuss the jet function. We define 

(0|T{(7-ec,.c)«(x)(e.,.e7-)/3(0)}|0> = ^ J ^e"^^V^(A:) , (117) 

where a and (3 are Dirac indexes. So far we have not specified the quantization frame. If 
we quantize in the equal-time frame we obtain at tree level 

J- = 2j- ^-^—=-y~J-. (118) 

If we quantize in the light-front frame with our standard gauge fixing prescription = 0, 

2 2 

the ^c,hc field becomes a constraint (the factor ^ ~ ^ is treated as a correction when 
quantizing for consistency). Either way, the imaginary term (the one that appears in the 
decay rate) reads (to be kept in mind that k~ ~ itlq > 0): 

ImJ- = -27i-f-S{k~^ -ml/k-) . (119) 

The soft function correlator reads 

HQ\Qlr{x-)P[exp{ig j^^ dz- A+{z-))]QM\H<^ ^2 j dp^e^^^^' S- (p^) , (120) 

where S~ is usually named the shape function. 

We are now in position of writing dT^~ydx in terms of . We obtain 



dr(-) 1 1 



dx Mhq 2(27r)x 27r 

X J d^xe-^'^e'""^'''' J dp-^e'"-^S-{p^) j (^^e-^'^'^mJ- . (121) 
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S{x) 




0.2 0.4 0.6 0.8 1 



Figure 12: Plot of the shape function with x = p'^ / P^^ for the values rris = (3 and 
rriQ = 10/5 (solid line) and nis = (3 and mq = 3(3 (dashed line). Strictly speaking the 
shape function is singular for x ^ 0. Nevertheless, this can only be seen by the eye for 
low values of niq. 



By doing the x'^ integration and using the fact that the x'^ variations of the field are 
small in comparison with mq (already used in the above expression), one can set k~ = mq 
in ImJ~ , which is correct at leading order in the 1/mq expansion. 

At this stage there is an strong simplification if we work in the light-front presented 
in sec. El with gauge fixing A'^ = 0, since this allows us to write the shape function in 
terms of the Hq-meson wave function squared, which now reads (p"*" + mqv^ = xM^q ) 

^ Ml^ x^ -M^/^^^'^^- ^^^^^ 

Note that we can do that because the shape function is computed at "equal" times (i.e. 
for "x"*" = constant") and we can use the free expressions for the heavy quark fields. In 
the last equality we have used the fact that terms of order 1 — x go like {3 /mq (at least 
when we compute moments). We show how the shape function looks like in Fig. ^| It is 
remarkable how similar it is to what one would expect in four dimensions. 

Actually one could do a similar analysis in four dimensions as far as one neglects 
higher Fock components of the Hq meson state. In this approximation one could relate 
the shape function with the square of the Hq wave function and see what is the impact, 
for instance, in the analysis of Ref. |25j . 

We could also play the same game for dT'^^^dx. We should then redo the construction 
of the effective theory, since the hard-coUinear quark would actually go in the opposite 
direction. Everything would work analogously changing J~ and S~ by and . 
reads equal to J~ changing 7" 7"*" and ^ k~ . Therefore ImJ^ and ImJ~ produce 
the same delta function (as far as the x variable is concerned). The definition of would 
come from Eq. ()120j) . changing 7" — 7"*" and x^ ^ x~ . This has consequences when 
working in the light-front frame. If we work in the gauge = everything is completely 
analogous to the computation of dV^^ydx, changing + ^ — everywhere. In particular 
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this implies that = S^. It is more interesting to consider our standard gauge fixing 
= 0. In this situation the hard-collinear field is not a constraint anymore, but the 
price to pay is that cannot be easily computed because the fields act at different times. 
Therefore, it is welcome that we can obtain 5*"^ from parity arguments. 
Our final results for ImT and 4^ read 

ax 

I^y(..eeZe.eO ^ ^^^^2^^ ( O; + ^ ) / , (123) 




j-p {tree level) -r -i ^2 

^ = ^—^(MH^x?2lmT^''''^''''^'^ (124) 



^ [x^ 



where in both expressions, we have included the subleading kinematical corrections, and 
those due to the mass of the hard-collinear, which are also parametrically subleading. 
The moments then read 

^,j.„.,„,„ _ I (X). (125) 

We first note that this expression agrees with the result shown in Eq. ()83p at LO in I/itlq. 
The 0{rri'l/rn?Q) corrections are also correctly incorporated. The first discrepancies are of 
0(/3^/mQ). We need a one-loop analysis to incorporate them, which we postpone to the 
next section. 

Finally, we would like to emphasize that the effective field theory computation has lost 
the information of the final bound state. This can be visualized by comparing Eq. 
with the exact hadronic expressions, which consist of a discrete sum over resonances. 



5 Soft-(Collinear) Effective Theory \YiD = l + l 

In the previous section we have worked out SCETI and the differential decay rate at 
leading order in A. The procedure seems to be quite cumbersome. In principle, it could 
be possible to extend the previous analysis to higher orders in A. Nevertheless, the process 
becomes tedious and, in principle, other modes should be included in the theory. This 
approach does not explicitly profit from the fact that hard-collinear modes appear from 
a very specific interaction (the weak one) and that the interaction takes place at very 
short "times" in the x'^ axis. Instead we will follow here an alternative approach and 
we will derive the effective theory by integrating out any (light) particle with P~ ~ uiq. 
This will make the effective field theory non-hermitian, introducing imaginary terms in 
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Figure 13: Symbolic plot to represent the matching from QCD to the effective theory of the 
(tree-level) diagram, the imaginary part of which produces the decay of the heavy quark to 
a hard-collinear quark and the particle. The RHS of the figure represents the effective 
vertex in the effective theory. The shape of the effective vertex represents that it is local 
in but not in x^ . 



the Lagrangian. This should not be considered unusual since this also happens for non- 
relativistic effective theories. The idea, which we have already mentioned previously, is 
that diagrams for which the P~ momentum flow are of the order of mq are local with 
respect the x^ axis in the light-cone frame and, once we have chosen the quantization 
frame to be = constant, we can construct the corresponding effective Hamiltonian. 
Therefore, the effects due to the hard-collinear fields can effectively be reproduced in the 
effective theory by a local vertex in x^ although not local in the other components (see 
Fig. fTT?|l . We should stress that at the end of the day we only want to describe the 
differential decay rate at leading order in the weak interactions. Therefore, they can be 
described by a purely imaginary vertex interaction. 

The effective degrees of freedom of our theory will be only soft ones (at least to the 
order to which we will work here). This considerably simplifies our approach compared 
with others where the derivation of the effective theory is performed investigating all 
possible modes existing in the theory. Actually one should not call this effective theory 
" soft-collinear" , since the collinear modes do not appear in the theory anymore, rather 
one may call it HQET but with some non-hermitian terms. 

The fact that our interactions will be almost local in a;"'"-times will also have important 
consequences in how the computation is performed. In principle, once all the dependence 
in has disappeared from the interaction^^, the interaction only depends on Pq and 
g"*". This means that one can use the free-field expressions for the fields in the effective 
interaction. Note that this is so even if the momentum P^ flowing in the diagram is 
small. We also have to take into account that we know the explicit expressions of the 
bound states in terms of the free-field expressions of the field. Therefore, we will be able 

-"^^More precisely, one can replace and by a function of niQ, Pq and g+ using the equations 
of motion up to residual 2;+ "time-dependent" terms, which can be expanded in an expansion in I/toq. 
One could get rid of these time-dependent terms systematically using field redefinitions. Nevertheless, in 
this paper, we do not reach enough precision to worry about this problem. 
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Figure 14: Tree level diagram, the imaginary part of which gives the leading contribution 
to the effective vertex in the effective theory, Eq. il29\) . 



to obtain explicit expressions of the matrix elements in terms of the 't Hooft wave function 
of the bound states, once we project the effective interacion to the physical states. It is 
then when the non-perturbative dynamics appears. 

An important point here is that the effective vertex can be obtained using perturbation 
theory, order by order in Since jS"^ has dimensions in two dimensions, the inclusion of 

terms produces subleading effects in the OPE. This is so because, even if ~ /5, Pj^ 
typically appears in the combination (except for the leading term) jS'^/P^, which could be 
interpreted as ~ Therefore, if we restrict ourselves to the kinematical situation 

Px 3> the effective interaction can be obtained using perturbation theory. This is 
what we will do in what follows. 

Our starting point is the Lagrangian of QCD where the gluons and the ■?/'_ component 
have been integrated out, i.e. Eq. (jH), coupled to the electroweak effective vertex pTj) . 
We then want to integrate out any (light) degree of freedom with P~ ~ mg and build 
an effective theory with only soft (light and heavy) quarks. Then we match the effective 
theory onto QCD. 

At tree level we only have to compute one diagram, which we have to match onto the 
effective vertex (note that throughout we only demand the imaginary piece to be equal). 
This is simbolically displayed in Fig. 

So far we have not specified neither the gauge nor the quantization frame. For the 
gauge fixing we follow our standard prescription of = 0. Following the notation of 
the previous section, at lowest order, the intermediate hard-coUinear field becomes a free 
particle and its propagator reads (A;~„ = m^/fc"*", k'^j^ = k^) 

T- - hn+m _ - f l^k-^ 7-fc+ \ _ 1 

^rn^ 1 _ - 7- 

k~ 

if we quantize in the light-front frame, and 



^ + m 

m? + ie ' ' k^ — 



J- = 7- .. ^T, .X = 27-7^^ , (127) 
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Figure 15: One loop diagrams (plus their symmetries), the imaginary part of which con- 
tribute to the effective vertex in Eq. iiy(J\) . 



if we quantize in the equal-time frame. We note that the the imaginary term (the one 
that appears in the decay rate) is equal in the light-front or equal-time quantization frame 
(to be kept in mind that k~ ~ tuq > 0): 



ImJ 



(128) 



Putting everything together, the contribution to the vertex reads (where we have used 



the free equation of motion k 



Q 



m 



effective vertex ~ Im 



Q 



(129) 



At 0(/3^), the following diagrams have to be considered (Fig. IT^ . Their net effect 
is to renormalize the masses of the hard collinear and heavy quark. The first diagram 
renormalizes the hard-coUinear mass that appears in the hard-coUinear propagator. The 
effect of the second diagram is to renormalize the masses of the vertex. To simplify the 
expression we keep subleading terms in the l/mq expansion. 

By adding all the terms, the effective interaction can be written as {P^ = Pq — q^) 



effective vertex 



Im 



^Ir ^q,r 



p2 



-VT- T—rO 



m 



(130) 



(P'c. 



P 



m: 



Q,R, 



where with the precision of the calculation, we have replaced mg by mg ^ within the 
delta. 

We are now in the position to write the effective Lagrangian, which, in fact, is the 
HQET Lagrangian adding the effective vertex. 



C- = Chqet + Im[Cj] 



[131) 
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and the effective vertex reads 




In order to Iceep the expression for the effective vertex more compact, we have written it 
in terms of the field had we written it in terms of Q+^i there would be a shift in the 
derivatives (for instance id^ rnqv^ + id^). We note that the Lagrangian is local in . 



5.1 Semileptonic differential decay rate 



At this stage we can compute the semileptonic decay. Actually we will have to compute 
the imaginary part produced by the effective vertex 



ImT, 
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The differential decay rate then reads 

1 1 



dx 



Mhq 2{2tt)x 2tt 



{MHQxf2lmT,ff 



(134) 
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Let us note that we can also rewrite this result (actually one equivalent with the precision 
of our calculation) in terms of shape and jet functions as in Eq. ()12H) . For the shape 
function, Eq. ()122p . we have to replace mg — > mg /j and for the jet function, Eq. pi8|) . 
we have to replace rric mc,R- 

The expressions for the moments read 
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'PHc 
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This expression allows us to compare the hadronic and the effective field theory compu- 
tations. It is quite remarkable that this result agrees with the hadronic result obtained in 
Eq. (jHHj) . Therefore, we also obtain the same expressions for M^^^, Eq. (jHll), valid in the 
OPE kinematic region, and for M^^^'^^ , Eq. (j85|) . valid in the SCETI kinematic region. 
The comments about the precision of the expressions made there also apply here. The 
expression for M^^^ is correct up to, and including, ©(l/m^) corrections. The expres- 
sion for Mfp^'^^ is correct up to, and including, 0{l/mQ) corrections. The conclusion is 
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that, if working with moments, we do not see duahty violations with the precision of our 
calculation. 

It is interesting to compare with the computation of the moments made by Bigi et 
al. For the total decay width. Mi, we obtain exactly the same analytic expression. 
Actually, this result is not affected by the radiative corrections (we obtain the same 
result at tree level or one-loop). In our opinion, this explains the agreement, since the 
computation of Bigi et al. has (effectively) been done at tree level. For Mq^2, we need 
to perform the one-loop computation in order to get agreement with their results with 
©(l/mg) precision. This is quite remarkable, since in the kinematical situation they 
choose, q'^ = 0, it is argued that there is a non-renormalization theorem for the current, 
whereas in the kinematical situation q~ = this is not true and the one-loop corrections 
to the vertex shown in Fig. ^Jhave to be included. 

We are now in the position to perform a numerical analysis of the corresponding 
expressions obtained from effective field theories with perturbative factorization and to 
compare them with the hadronic ones. We first consider the differential decay rate. As we 
have already discussed throughout the paper, the direct comparison between the partonic 
and hadronic result is not possible since the first is an smooth function in x, whereas the 
second consists of a sum over deltas. Nevertheless, the layer-function approximation gives 
us a qualitative way to compare the hadronic matrix elements with the computation using 
the effective theory, since it allows us to write the hadronic matrix elements in terms of 
only the Hq meson wave-function. Naively, one could make the following assignment to 
the hadronic matrix element from the effective field theory computation 



where the normalization has been adjusted to agree with the layer-function result. The 
LO result agrees with the LO layer-function expression, Eq. (jfifjj) . Nevertheless, this is 
not so for the subleading one. We will elaborate on this issue at the end of this section. 
Although numerically the effect is not very important (we would obtain very similar plots 
to those obtained in Figs. [7|and|Hl), it is important from the conceptual point of view. 

We now move to the comparison of the moments. In this case, the comparison is more 
sound. For all the numerical checks we have set G^/(47r) = 1. We have compared the 
perturbative and hadronic result. The agreement is very good up to very high moments. 
We show this comparison in Fig. I16r^. We also consider how convergent is the expansion 
either if we work in the strict OPE or SCETI regime. We show the convergence of the 
OPE expansion in Fig. El For the OPE limit we can see that the breakdown of the 

^''We also make the comparison for the set of masses used in Ref. to iUustrate how weU our 
expressions would work for an hypothetical b — s- civ decay. 




(136) 
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Figure 16: Difference between the hadronic, Eq. ^7^ , and perturbative, Eq. result 
for the moments (divided by the hadronic result). The first figure is for the values mq = 
10(3, rric = (3 and rUg = (3, the second is for the values mq = 15/5, nic = 10/5 and 
rus = 0.56/5 and the third for the values mq = 3/5, nic = P and rUg = (3. 



agreement with the hadronic result appears earher. This is to be expected since there is 
a new scale N/mq, which is not resummed. In any case the precision is better than 5% 
for N below 6 at NNLO. We also perform the analysis in the SCETI region. We show 
the plot in Fig. ^1 As expected the convergence improves over the OPE evaluations 
and they are optimal for values of ~ rnq/P as expected. For very large values of A^, 
they blow up as expected. We also consider the dependence on the heavy quark mass of 
these results. For that we repeat the same analysis for mq = 3, which is somewhat the 
limiting case of validity of our results. Overall, we get a similar picture than before but 
with worse convergence (actually in the OPE region the results are barely convergent) 
and the perturbative results are only reliable for lower moments, again as expected. The 
effect of the one-loop corrections (the renormalization of the masses) is small. 

'■ — 2 — I improves the numerical agreement with the 

^Q,R J 

hadronic expression but we should remind that we cannot claim better accuracy than 
^ 1 — A^ I'"" . This effect is particularly important if the mass of the hard-coUinear is 

Q,R 

large. 
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Figure 17: Difference between the hadronic, Eq. (|7ij| ), and perturbative, Eq. \84^ , result 
for the moments (divided by the hadronic result) in the OPE limit. The dashed line for 
the LO result, the dash-dotted line for the NLO result, and the solid line for the NNLO 
result. We take the values mq = lOP, rric = P and nis = P for the first figure and 
ruQ = SP, rric = P and rus = P for the second figure. We use the values {t) = 1.73/5 and 
{t^) = 3.99/5^, which can be checked with the sum rules of Ref W^ . 



Rel . error Rel . error 




N 



Figure 18: Difference between the hadronic, Eq. and perturbative, Eq. \85]) . result 

for the moments (divided by the hadronic result) in the SCETI limit. The dashed line 
for the LO result, and the solid line for the NLO result. We take the values mq = lOP, 
rric = P and = P for the first figure and rriq = SP, = P and = P for the second 
figure. 
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We now consider THgiy)- We compare the hadronic result versus the prediction from 
effective field theories in Fig. ^1 We can see a good quantitative agreement between both 
lines but it is not possible to perform a point-to-point comparison because the hadronic 
result is a step-function whereas the perturbative computation is an smooth function in 
y. Therefore, it is difficult to quantify the error, which typically is of the order of the 
difference between the n and n + 1 state contribution to Th^ (y) . 

To make more evident this problem, we now consider the average 



x„+Sx 

—dx , 



Xn — Sx 



dx 



(137) 



where 6x is bounded to be small enough that only one resonance contributes to the 
integral. One could believe that the duality violations can be smoothed in this way. This 
definition provides the closest possible thing to a point-to-point comparison between the 
perturbative and hadronic result. Since we are now able to perform analytic computations 
both for the hadronic and for perturbative result, we are in the position to quantify this 
statement (we restrict to the kinematic regime where ^ (3'^). We obtain the following 
result from the hadronic computation 



-Sx dx 



dx 



From the perturbative computation we obtain 
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(139) 

We can see that both expressions are different, even at leading order, due to the normaliza- 
tion. Only if we fine tune 5x to a very specific value (actually, the one chosen in Ref. [S]), 
we get agreement between both expressions, and, even then, the subleading corrections 
are different^^. In order to fine tune the value of 5x, we should have a good knowledge 
of the non-perturbative spectrum (this is certainly difficult in the four- dimension case, 
although one can always assume a linear regge behavior), which is something that we 
do not expect can be achieved from perturbation theory. On the other hand, as far as 
5x is independent of n, one obtains this equality at leading order (subleading corrections 
remain to be different) 



Jxn+i — 5x dxlhadr. 
1 



X„+Sx I 

x„—5x dx I hadr. 



dx 



p+^+ffl dx 

•> Xn+i—ox dx \pert. 

p+f^l dx 

Jxn—ox dx Ipert. 



(140) 



without fine-tuning 6x. 



^^Actually, these small differences for the subleading corrections are crucial to get agreement for the 
moments at 0(1/toq). 
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Figure 19: Plot of T Hqiv) using the hadronic or perturbative expression of dV/dx. The 
smooth curve represents the perturbative prediction and the other line is the hadronic one. 
We take the values niq = 10(3, rric = (3 and = (3 for the first figure and mq = 3(3, 
rric = (3 and ntg = (3 for the second figure. 

6 Conclusions 

We have studied QCD in 1+1 dimensions in the large Nc hmit using hght-front Hamilto- 
nian perturbation theory in the 1/A^c expansion. We have used the formahsm developed 
to exactly compute hadronic transition matrix elements for arbitrary currents at lead- 
ing order in l/Nc- We have compared with previous results found in the literature. We 
have computed in two alternative ways the semileptonic differential decay rate of a heavy 
meson and its moments using the previously computed hadronic matrix elements. They 
yield very different expressions for the differential decay rate, which should be equal by 
parity invariance. This has lead to the derivation of non-trivial equalities between matrix 
elements, which we have checked numerically. Some partial analytic checks have also been 
done. 

We have then focused on the kinematic regime where the OPE (A^ ~ x ~ 1) or 
SCETI (1 — a; ~ 1/A^ ~ P/'^q) can be applied. This means n S> 1, where n is the 
principal quantum number of the final hadronic state. This has allowed us to use the 
properties of the final hadronic bound state in those kinematical regimes by using the 
layer function |7pi8j . The resulting expressions are suitable to a more direct connection 
with the computation obtained using perturbative factorization. We have then obtained 
expressions for the moments using the EulerMcLaurin formula, within an expansion in 
l/mg, either in the OPE or in the SCETI regions. In the first case up to order l/rriQ and 
in the second up to order l/mg. We have also checked that these results agree with the 
expressions of the sum rules for = 0, 1,2 obtained in Ref. |1] up to order l/rriQ (we 
have also checked them by direct computation). 

We have also studied the differential decay rate using effective theories with perturba- 
tive factorization. We have first derived SCETI at leading order in 1+1 dimensions, and 
applied it to the differential decay rate. We have seen that there is an strong simplification 
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working in the light-cone quantization frame, which allows us to relate the shape function 
with the wave-function of the bound state (in the large N^.). We have then dwelt further 
on the issue of finding the optimal effective theory to describe the differential decay rate 
in the OPE and SCETI kinematical regime. We see that, at least in two dimensions, it 
appears to be more efficient to integrate out the hard-coUinear modes and only keep soft 
degrees of freedom. This takes advantage of the fact that the hard-coUinear interactions 
only appear in a very specific way, emanating from the weak interaction, whereas in the 
standard construction of SCETI one works out all possible couplings. The resulting effec- 
tive theory is equal to HQET plus some imaginary terms, which describe the differential 
decay rate. This effective theory is more efficiently implemented in the hght-front quanti- 
zation frame, where it becomes "local" in the "x"*"" quantization frame, and it is suitable 
for computations in a Hamiltonian formulation. We have then obtained the differential 
decay rate at one-loop. We would also like to remark that some of these ideas could 
also apply to the four-dimensional case, where it also seems possible to obtain a "local" 
effective interaction in the "x"*"" quantization frame. 

At the end of the day, we have been able to obtain expressions for the moments using 
effective field theories with perturbative factorization with relative accuracy of 0{f3'^/mQ) 
in the kinematical regime where the OPE can be applied, and with relative accuracy of 
0{j3/mQ) in the kinematical regime where non-local effective field theories can be applied. 
These expressions agree, within this precision, with those obtained from the hadronic ex- 
pressions using the layer-function approximation plus Euler-McLaurin expansion (taking 
into account the numerical agreement for the radiative correction). 

Numerically very good agreement for the moments between the exact result and the 
result using effective field theories with perturbative factorization is obtained. For the dif- 
ferential decay rate it is also possible to perform a sort of comparison between the hadronic 
and perturbative result based on the layer-function approximation. The agreement is very 
good if we do not approach too much the limit x — > 1. 

There is still the issue of the theoretical error. As we have mentioned throughout 
the paper, it is not possible to make a quantitative comparison (with errors) between 
the hadronic and perturbative differential decay rate, since one is represented by a sum 
of deltas whereas the other is a smooth function in x. Therefore, effective field theories 
with perturbative factorization are not suitable to predict the differential decay rate on a 
point-to-point basis. Note that this comment also applies to THqiv) and, in principle, to 
other arbitrarily smeared functions. This is best illustrated in the derivation of SCETI 
at leading order where, after field redefinitions, one is lead to an effective theory where 
the hard-colhnear is a free field. Therefore, it can never build a bound state, which is 
what is observed in the hadronic differential decay rate. Effective field theories with 
perturbative factorization can only hopefully be a good aproximation (or at least a good 
starting point) for inclusive observables on which one averages over a large fraction of the 
final bound states (this is a handicap from the experimental point of view, since one has 
to know the differential decay rate for arbitrarily large momenta). This motivates the 
use of moments^^ for the comparison, since they may lead to a "quantitative" comparison 

^^When we perform the comparison among moments the difference is "less severe" than using directly 
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between hadronic and perturbative results. With the precision obtained in this paper 
we obtain a perfect match between the hadronic and perturbative results. Nevertheless, 
to have a more rigourous handle of the errors, one should be able to quantify the error 
produced by using the layer-function (i.e. to consider subleading effects in the WKB 
approximation), as well as by using the Euler-McLaurin formula. Actually it is this 
last formula that allows to make quantitative the comparison of perturbative and non- 
perturbative results. We expect to come back to these issues in the future in order to try 
to find duality violations in the computation of the moments in the SCETI kinematical 
region. 
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A Conventions and notation 

In this appendix the conventions and notation that we use are presented. We define two 
light-like vectors (with the metric g^" = g ^ = 2 and zero elsewhere), 

n^: = (l,l), < = (!,-!), (141) 

where light-cone coordinates are defined in the usual way, 

x"*" = ■ X = + x^) , x~ = n_ ■ X = — x^) , (142) 

which imply that 



x° = - (x+ + a;~) , x^ = -(a;+-x), (143) 



and 



d^ = 2j^ = /^ + /^ = do+ds ~ , a+ = 2-^ = = do-d, ~ p+ , (144) 

OX+ ox^ ox-^ ox~ ox" ox-^ 



P+x- P'x+ 
P ■x = + , (145) 



d'^x = ^dx+dx-d^-^xs_ . (146) 



the differential decay rate. The point is to quantify the "less severe" at the parametrical level (although 
it can be done at the numerical level). For the moments iV = 0, 1, 2, this has already been done in Ref. 

El 
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For the Dirac algebra is useful to define the corresponding light-cone matrices 

7^+ = 7+, i_=r ■ (147) 

To have explicit expressions, it is useful to work with an explicit representation of the 
Dirac algebra. We will use the following Weyl-like representation for the Dirac algebra 

-'°-(! o ) ^'-(° o)' ''''' 

so that the corresponding light cone matrices are given by 

This way one can see the explicit effect of the projection operators (75 = 7''7^) 
_ 1 + 75 7V 1 . . 1 _ + /I 0\ 

_ 1-75 7°7" 1 . . 1 + - / 
^- = -^ = ^ = 4^+^- = r ^ =[0 1 

in the sense that if one splits a quark doublet on its two components 

^=(f,:)^ (151) 



then the projection operators act as expected 
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